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1.  Introduction 


This  report  presents  a  set  of  functions,  written  in  C++,  that  can  be  used  to  calculate  atmospheric 
conditions  (temperature,  pressure,  air  density,  and  speed  of  sound),  as  well  as  gravitational-field 
strength.  The  functions  are  based  on  the  atmospheric  model  presented  in  U.S.  Standard 
Atmosphere,  1976 } 

A  summary  sheet  is  provided  at  the  end  of  this  report.  It  presents  the  yAtmosphere  namespace, 
which  contains  the  five  functions  that  are  described  in  detail  in  this  report.  Also  presented  is  an 
example  that  can  be  used  to  recreate  portions  of  the  “Main  Tables”  presented  in  U.S.  Standard 
Atmosphere,  1976. 


2.  Summary  of  Constants 


The  following  is  a  list  of  constants  that  are  used  by  the  code  presented  in  this  report.  All  are 
from  U.S.  Standard  Atmosphere,  1976.  Page  references  refer  to  the  source  document. 


R*  =  83 14.32  N  m/(kmol  K) 

universal  gas  constant 

(page  3) 

r0  =  6,356,766  m 

effective  radius  of  the  earth 

(page  4) 

'N 

II 

O 

o 

specific  heat  ratio 

(page  4) 

Tq  =288.15  K 

standard  temperature  at  sea  level 

(page  4) 

g'Q=  9.80665  m2/(s2m') 

geopotential  constant 

(page  8) 

gQ  =9.80665  m/s2 

gravitational-field  strength  at  sea  level 

(page  8) 

r  —  1  mVm 

unit-converting  constant 

(page  8) 

M0  =  28.9644  k^kmol 

mean  molecular  weight  of  air 

(page  9) 

P0  =  101325.0  Pa 

standard  pressure  at  sea  level 

(page  12) 

1  National  Aeronautics  and  Space  Administration.  U.S.  Standard  Atmosphere,  1976\  NASA-TM-X-74335;  U.S.  Government 
Printing  Office:  Washington,  DC,  October  1976. 
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3.  Summary  of  Equations 


The  following  is  a  list  of  equations  that  are  used  by  the  code  presented  in  this  report.  All  are 
from  U.S.  Standard  Atmosphere,  1976.  Equation  numbers  refer  to  the  source  document. 


8  =  8  c 


r0+z 


gravitational-field  strength 


(equation  17) 


rn  +  Z 


geopotential  height 


(equation  18) 


TM=TM,b+LM'b(H-Hb) 


molecular-scale  temperature 


(equation  23) 


Tmj,  +!«(»-«,) 


-So  M0(H-H„) 


for  LM.b  =  0 


for  lm.„  *  0 


pressure 


(equation  33ab) 


air  density 


(equation  42) 


yR  tm 


speed  of  sound 


(equation  50) 


4.  Augmented  Table  4 


Along  with  equations  23  and  33ab,  U.S.  Standard  Atmosphere,  1976  specifies  the  use  of 
tabulated  values  to  calculate  temperature  and  pressure.  Values  for  Hh  and  LM  b  are  listed  in  the 

source  document’s  table  4.  Values  for  TM  b  and  Ph  can  be  found  by  treating  equations  23  and 
33ab  as  recursive,  with  T0  and  P0  from  section  2  used  as  initial  values.  In  practice,  it  is 
computationally  efficient  to  precalculate  values  for  TM  h  and  Ph .  Table  1  lists  the  values  for  Hh 
and  LM  b  that  were  originally  listed  in  the  source  document’s  table  4,  as  well  as  precalculated 
values  for  TMh  and  Pb. 
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Table  1.  Tabulated  values  for  H h  ,  LM  b  ,  TM  h ,  and  Pb  . 


b 

Hb 

(m1) 

P  M,b 

(K/m’) 

T M,b 

(K) 

Pb 

(Pa) 

0 

0 

-0.0065 

288.15 

1.01 325000000000E+05 

1 

11000 

0 

216.65 

2. 263206397 34629E+04 

2 

20000 

0.001 

216.65 

5. 47488866967777E403 

3 

32000 

0.0028 

228.65 

8. 680 1 86847 55228E+02 

4 

47000 

0 

270.65 

1. 10906305554966E+02 

5 

51000 

-0.0028 

270.65 

6. 6938873 1186873E+01 

6 

71000 

-0.002 

214.65 

3.95642042804073E-fO0 

7 

84852 

0 

186.946 

3.73383589976215E-01 

5.  Calculating  Temperature  and  Pressure 


The  following  method  is  valid  for  -  5,000  m  <  z  <  86,000  m .  Begin  by  calculating  geopotential 
height  ( H )  as  a  function  of  altitude  ( z  )  using  equation  18.  Next,  find  Z?  I  H ^  ^  H  Hb+l  ■  If 
H  <  0  m' ,  use  b  =  0 .  If  H  >  8,4852m' ,  use  b  =  7  . 

To  find  molecular- scale  temperature,  substitute  values  for  Hh,  LM  h ,  and  TM  b  into  equation  23. 
To  find  pressure,  substitute  values  for  Hb ,  LMh,  TM  b ,  and  Ph  into  equation  33ab. 


6.  C++  Implementation 


6.1  TABLE4  Array 

The  values  for  H h ,  LM  h .  TM  b ,  and  Pb  are  stored  in  the  TABLE4  array. 

TABLE4  Code 


const  double  TABLE4[8][4]={//< 

- TRANSITI0N  points 

00000 

3 

-0.0065 

3 

288.150 

3 

1 . 01325000000000E+5 

,// 

FOR  PRESSURE  & 

11000 

» 

0 . 0000 

3 

216.650 

3 

2 . 26320639734629E+4 

,n 

TEMPERATURE  VS 

20000 

> 

0 . 0010 

3 

216.650 

3 

5.47488866967777E+3 

,// 

GE0P0TENTIAL 

32000 

3 

0.0028 

3 

228.650 

3 

8 . 68018684755228E+2 

,// 

ALTITUDE  CURVES 

47000 

3 

0 . 0000 

3 

270.650 

3 

1 . 10906305 5 54966E+2 

>n 

[table  4] 

51000 

3 

-0.0028 

3 

270.650 

3 

6 . 69388731186873 E+l 

,// 

(3RD  COLUMN  IS 

71000 

3 

-0 . 0020 

3 

214.650 

3 

3 . 95642042804073E+0 

,n 

TEMPERATURE, 

84852 

3 

0 . 0000 

3 

186.946 

3 

3 . 73383 58997621 5 E  -1 

n 

4TH,  PRESSURE) 

};// - YAGENAUT@GMAIL 

.  C  01^1°^  '■v 

_AST  ~UPDATED~15MAR2013~ - - 

3 


6.2  TemperatureO  Function 


The  Temperature()  function  uses  the  method  described  in  section  5  to  calculate  molecular-scale 
temperature  as  a  function  of  altitude.  For  altitudes  up  to  80,000  m,  molecular- sc  ale  temperature 
is  equal  to  absolute  temperature.  At  80,000  m,  molecular-scale  temperatures  begin  to  differ  from 
absolute  temperatures.  At  86,000  m  (the  limit  of  this  model),  U.S.  Standard  Atmosphere,  1976 
states  that  the  difference  between  molecular-scale  temperature  and  absolute  temperature  is 
0.0787  K. 

TemperatureO  Code 


inline  double  Temperature(//<==================================TEMPERATURE  (K) 

double  z){//< . ALTITUDE  (m)  (T  IS  VALID  FOR  -5,000  m  <  z  <  86,000  m) 

double  H=z* 63 56766/ ( z+63 56766 ) ; / / . [equation  18] 

int  b;/*<-*/for(b=0; b<7;++b)if (H<TABLE4[b+l] [0] ) break; 

return  TABLE4[b] [2]+TABLE4[b] [l]*(H-TABLE4[b][0]);// . [equation  23] 

}// —  ~~YAGENAUT  (5)GMAI  L .  L  AS  T  ~UPDATED~1 5  MAR  2  0 1 


TemperatureO  Parameters 

z  z  is  a  height  above  sea  level  (in  meters).  Output  of  the  TemperatureO  function  is  valid 
for  -5000  m  <  z  <  86,000  m.  Output  for  z  =  -6356766  m  (i.e.,  the  center  of  the  earth) 
is  undefined. 

TemperatureO  Return  Value 

The  return  value  of  the  TemperatureO  function  is  the  molecular-scale  temperature,  in  Kelvin,  at 
the  given  altitude  (z).  See  the  beginning  of  this  section  for  a  note  on  the  minor  difference 
between  molecular- scale  temperature  and  absolute  temperature. 

6.3  PressureO  Function 

The  Pressure()  function  uses  the  method  described  in  section  5  to  calculate  atmospheric  pressure 
as  a  function  of  altitude. 

PressureO  Code 


inline  double  Pressure(//<=======================================PRESSURE  (Pa) 

double  z){//< . -ALTITUDE  (m)  (P  IS  VALID  FOR  -5,000  m  <  z  <  86,000  m) 

double  H=z*6356766/ (z+6356766) ; / / . [equation  18] 

int  b;/*<-*/for(b=0; b<7; ++b) if (H<TABLE4[b+l] [0] ) break; 

double  C=- .0341631947363104;// . C  =  -G0*M0/RSTAR  [pages  8,9,3] 

double  Hb=TABLE4[b] [0] , Lb=TABLE4[b] [1] ,Tb=TABLE4[b] [2] , Pb=TABLE4[b] [3]; 
return  Pb*(fabs(Lb)>lE-12?pow(l+Lb/Tb*(EI-Hb) ,C/Lb)  :exp(C*(EI-Elb)/Tb)); 

}//  Y  AG  E  N  AU  T  (5)GMA  I L .  C  LAST  ~U  P  DAT  E  D~  1 5  MA  R  2  0 1 
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PressureO  Parameters 

z  z  is  a  height  above  sea  level  (in  meters).  Output  of  the  PressureO  function  is  valid  for 
-5000  m  <  z  <  86,000  m.  Output  for  z  =  -6356766  m  (i.e.,  the  center  of  the  earth)  is 
undefined. 

PressureO  Return  Value 

The  return  value  of  the  PressureO  function  is  the  pressure,  in  Pascals,  at  the  given  altitude  (z). 

6.4  Density ()  Function 

The  DensityO  function  uses  equation  42  to  calculate  air  density  as  a  function  of  temperature  and 

pressure. 

DensityO  Code 

inline  double  Density(//<=====================================DENSITY  (kg/mA3) 

double  !,//< . TEMPERATURE  (K)  (CALCULATE  T  USING  Temperature( ) ) 

double  P){//< . -PRESSURE  (Pa)  (CALCULATE  P  USING  Pressure^) 

return  P* .00348367635597379/T; // . [equation  42] 

J I  j  y  AG  E  N  AUT  |S)GMAI  | t  G  LAST  ^  LJ  P  DAT  E  D*^  1 5  MA  R  2  0 1 3 


DensityO  Parameters 

T  T  is  a  local  temperature  (in  Kelvin).  Values  for  T  can  be  determined  using  the 
Temperature()  function.  Output  for  T  =  0  K  is  undefined. 

P  P  is  a  local  pressure  (in  Pascals).  Values  for  P  can  be  determined  using  the  PressureO 
function. 

DensityO  Return  Value 

The  return  value  of  the  DensityO  function  is  the  local  density  of  air,  in  kilograms/cubic  meter. 

6.5  SpeedofSoundO  Function 

The  SpeedofSoundO  function  uses  equation  50  to  calculate  speed  of  sound  as  a  function  of 
temperature. 

SpeedofSoundO  Code 

inline  double  SpeedofSound(//<============================SPEED  OF  SOUND  (m/s) 

double  T) {//< . TEMPERATURE  (K)  (CALCULATE  T  USING  Temperature( ) ) 

return  sqrt(401 . 87430086589*T) ; // . [equation  50] 

}//  rsjrsjrsjrsj  YAGENAUT@GMAIL.COM - LAST-UPDATED-15MAR2013 - 


SpeedofSoundO  Parameters 

T  T  is  a  local  temperature  (in  Kelvin).  Values  for  T  can  be  determined  using  the 
Temperature()  function. 
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SpeedofSoundO  Return  Value 

The  return  value  of  the  SpeedofSoundO  function  is  the  local  speed  of  sound,  in  meters  per 
second. 

6.6  Gravity ()  Function 

The  GravityO  function  uses  equation  17  to  calculate  the  gravitational-field  strength  near  the 
surface  of  the  earth  as  a  function  of  altitude. 

GravityO  Code 

inline  double  Gravity(//<==================ACCELERATION  DUE  TO  GRAVITY  (m/sA2) 

double  z){//< . ALTITUDE  (m)  (g  IS  VALID  FOR  -5,000  m  <  z  <  86,000  m) 

return  9 . 80665*pow(l+z/6356766, -2) ; // . [equation  17] 

y  I  j  'nj/v/vpu  y  AG  E  N  AUT  (SGMA  XL*  C  0MruAjrwfvrw,v,v,Njrvrw,v,v,Nj,vrw,>j,v,Njrwrvrw,v,v,Njrw  LAST  ^  LJ  P  DAT  E  D*^  1 5  MA  R  2  0 1 3''"'"'"'"'"'' 


GravityO  Parameters 

z  z  is  a  height  above  sea  level  (in  meters).  The  output  of  the  GravityO  function  is  valid 
for  -5000  m  <  z  <  1,000,000  m.  Output  for  z  <  -6356766  m  is  undefined. 

GravityO  Return  Value 

The  return  value  of  the  GravityO  function  is  the  gravitational-field  strength,  in  meters  per  second 
squared. 


7.  Summary 


A  summary  sheet  is  provided  at  the  end  of  this  report.  It  presents  the  yAtmosphere  namespace, 
which  contains  the  five  functions  described  in  detail  in  section  6.  Also  presented  is  an  example 
that  can  be  used  to  recreate  portions  of  the  “Main  Tables”  presented  in  U.S.  Standard 
Atmosphere,  1976.1 

Some  values  that  are  calculated  using  the  example  code  differ  from  values  presented  in  U.S. 
Standard  Atmosphere,  1976  in  the  least-significant  figure.  For  example,  at  an  altitude  of  85,500 
m,  U.S.  Standard  Atmosphere,  1976  states  that  the  molecular-scale  temperature  will  be  187.920 
K,  while  the  example  code  predicts  that  the  molecular-scale  temperature  will  be 
187.919  K. 
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yAtmosphere  Summary 


y_atmosphere.h 


#ifndef  Y_ATMOS PH E R E_H_ 

#define  Y_ATMOSPHERE_H_ 

#include  <cmath>// . f abs  ( ) ,  exp( ) ,  pow( ) ,  sqrt (  ) 

namespace  yAtmosphere{//<======================================ATMOSPHERIC  MODEL 

/ /@@@@@@@@@@@@@@@  [BASED  ON  "U.S.  STANDARD  ATMOSPHERE,  1976"]  @@@@@@(S$@@@@@@@@ 
const  double  TABLE4[8] [4]={//<===============================TRANSITI0N  POINTS 


11000 

20000 

32000 

47000 

51000 

71000 

84852 


-0.0065 


0010 

0028 

0000 

0028 

0020 

0000 


0 


288.150 

216.650 

216.650 

228.650 

270.650 

270.650 

214.650 
186.946 


. 01325000000000E+5  ,// 
. 26320639734629E+4  , // 
5 . 47488866967777E+3  , // 
8 . 68018684755228E+2  , // 
1 . 10906305554966E+2  ,// 
6 . 69388731186873E+1  ,// 
3.95 64204280407 3 E+0  ,// 
3.73383589976215E-1  // 


FOR  PRESSURE  & 
TEMPERATURE  VS 
GEOPOTENTIAL 
ALTITUDE  CURVES 
[table  4] 
(3RD  COLUMN  IS 
TEMPERATURE, 
4TH,  PRESSURE) 


~LAST~UPDATED~15MAR2013~ 


};// - YAGENAUT@GMAIL.COM - 

inline  double  |T£Bi||©rattif’§(//<==================================TEMPERATURE  (K) 

double  z){//< . ALTITUDE  (m)  (T  IS  VALID  FOR  -5,000  m  <  z  <  86,000  m) 

double  H=z*6356766/(z+6356766) ; // . [equation  18] 

int  b;/*<-*/for(b=0;b<7;++b)if (H<TABLE4[b+l] [0] ) break; 

return  TABLE4[b] [2]+TABLE4[b] [l]*(H-TABLE4[b] [0] ); // . [equation  23] 

}// - YAGENAUT@GMAIL.COM - LAST-UPDATED-15MAR2013 - 

inline  double  pN?sSWPi(//<=======================================PRESSURE  (Pa) 

double  z ) {/ /< - ALTITUDE  (m)  (P  IS  VALID  FOR  -5,000  m  <  z  <  86,000  m) 

double  H=z*6356766/ (z+6356766) ; // . [equation  18] 

int  b;/*<-*/for(b=0;b<7;++b)if (H<TABLE4[b+l] [0] ) break; 

double  C=-. 0341631947363104;// . C  =  -G0*M0/RSTAR  [pages  8,9,3] 

double  Hb=TABLE4[b] [0], Lb=TABLE4[b] [1] ,Tb=TABLE4[b] [2] , Pb=TABLE4[b] [3] ; 
return  Pb*(fabs(Lb)>lE-12?pow(l+Lb/Tb*(H-Hb),C/Lb) :exp(C*(H-Hb)/Tb)); 

}// - YAGENAUT@GMAIL .  COM - LAST-UPDATED-15MAR2013 - 

inline  double  ftattsity(//<=====================================DENSITY  (kg/mA3) 

double  T, //< . TEMPERATURE  (K)  (CALCULATE  T  USING  Temperature( ) ) 

double  P){//< . PRESSURE  (Pa)  (CALCULATE  P  USING  PressureQ) 

return  P*  .00348367635597379/T;// . [equation  42] 

}// - YAG  E  NAUT@GMAIL.COM~ - LAS  T~U  PDAT  ED~15MAR2013 - 

inline  double  SpeedofSound(//<============================SPEED  OF  SOUND  (m/s) 

double  T){//< . . TEMPERATURE  (K)  (CALCULATE  T  USING  Temperature( ) ) 

return  sqrt (401. 87430086589*T);// . [equation  50] 

}// - YAGENAUT@GMAIL.COM - LAST-UPDATED-15MAR2013 - 

inline  double  l6paVitl(//<==================ACCELERATI0N  DUE  TO  GRAVITY  (m/sA2) 

double  z ) {/ /< - ALTITUDE  (m)  (g  IS  VALID  FOR  -5,000  m  <  z  <  86,000  m) 

return  9.80665*pow(l+z/6356766, -2);// . [equation  17] 

}// - YAGENAUT@GMAIL.COM - LAST-UPDATED-15MAR2013 - 


Example 


,  stdout);// . redirect  output 


#include  <cstdio>// . FILE, f reopen ( ), stdout, printf (),fclose( ) 

#include  "y_atmosphere.h"//Temperature(),Pressure(),Density(),SpeedofSound(), . . . 
int  main(){ 

using  namespace  yAtmosphere; 

FILE  *f=freopen("atmosphere.txt", "w1 
printf ( "#Z-Altitude  |  Temperature  | 

"|  Gravity\n"); 

printf ("#  (m)  |  (K)  | 

"|  (m/sA2)\n#"); 

for (int  i=0;i<79;++i) printf ("-"); 
for(int  i=0;i<1811;++i){ 
if ( i%10==0) printf ( " \n " ) ; 


I 


(Pa) 


Density 

(kg/mA3) 


|  Sound 
I  (m/s) 


double  z=-5000+i*50;/*->*/printf("%8.0f" ,z);// . altitude  (m) 

double  T=Temperature(z) ; /*->*/printf ( "%15 . 3f " ,T) ; // . temperature  (K) 

double  P=Pressure(z) ;/*->*/ printf ("%17.4E",P);// . pressure  (Pa) 

double  rho=Density(T, P) ; /*- >*/printf ( "%16 . 4E " , rho) ; // . density  (kg/mA3) 

double  Vs=SpeedofSound(T);/*->*/printf("%11.2f",Vs);// . sound  speed  (m/s) 

double  g=Gravity(z);/*->*/printf ( "%11.4f \n",g); }// . gravity  (m/sA2) 

fclose(f ) ; 
return  0; 

}// - YAGENAUT@GMAIL.COM - LAST-UPDATED-15MAR2013 - 


atmosphere 

.txt  (created  with 

the  above 

"Example”  code) 

#Z-Altitude  |  Temperature  | 

Pressure 

Density 

|  Sound 

Gravity 

#  (m)  | 

(K)  1 

(Pa) 

(kg/mA3) 

1  (m/s)  | 

(m/sA2) 

-5000 

320.676 

1.7776E+005 

1.9311E+000 

358.99 

9.8221 

-4950 

320.350 

1 . 7682E+005 

1 . 9228E+000 

358.80 

9.8219 

-4900 

320.025 

1. 7587E+005 

1 . 9145E+000 

358.62 

9.8218 

-4850 

319.699 

1.7493E+005 

1 . 9062E+000 

358.44 

9.8216 

-4800 

319.374 

1. 7400E+005 

1 . 8980E+000 

358.26 

9.8215 

-4750 

319.048 

1.7307E+005 

1 . 8898E+000 

358.07 

9.8213 

-4700 

318.723 

1. 7214E+005 

1 . 8816E+000 

357.89 

9.8212 

-4650 

318.397 

1.7122E+005 

1 . 8734E+000 

357.71 

9.8210 

-4600 

318.072 

1. 7030E+005 

1 . 8653E+000 

357.53 

9.8209 

-4550 

317.746 

1.6939E+005 

1 . 8572E+000 

357.34 

9.8207 

-500 

291.400 

1.0748E+005 

1 . 2849E+000 

342.21 

9.8082 

-450 

291.075 

1 . 0685E+005 

1 . 2788E+000 

342.02 

9.8080 

-400 

290.750 

1.0622E+005 

1 . 2727E+000 

341.83 

9.8079 

-350 

290.425 

1.0560E+005 

1 . 2667E+000 

341.63 

9.8077 

-300 

290.100 

1 . 0498E+005 

1 . 2607E+000 

341.44 

9.8076 

-250 

289.775 

1 . 0436E+005 

1 . 2547E+000 

341.25 

9.8074 

-200 

289.450 

1 . 0375E+005 

1 . 2487E+000 

341.06 

9.8073 

-150 

289.125 

1 . 0314E+005 

1 . 2427E+000 

340 . 87 

9.8071 

-100 

288.800 

1.0253E+005 

1.2368E+000 

340.68 

9.8070 

-50 

288.475 

1.0193E+005 

1 . 2309E+000 

340.49 

9.8068 

0 

288.150 

1.0133E+005 

1 . 2250E+000 

340 . 29 

9.8066 

50 

287.825 

1 . 0073E+005 

1 . 2191E+000 

340.10 

9.8065 

100 

287.500 

1 . 0013E+005 

1 . 2133E+000 

339.91 

9.8063 

150 

287.175 

9. 9536E+004 

1 . 2075E+000 

339.72 

9.8062 

200 

286.850 

9. 8945E+004 

1 . 2017E+000 

339.53 

9.8060 

250 

286.525 

9. 8358E+004 

1 . 1959E+000 

339.33 

9.8059 

300 

286.200 

9.7773E+004 

1 . 1901E+000 

339.14 

9.8057 

350 

285.875 

9.7191E+004 

1 . 1844E+000 

338.95 

9.8056 

400 

285.550 

9.6611E+004 

1 . 1786E+000 

338.76 

9.8054 

450 

285.225 

9. 6035E+004 

1 . 1729E+000 

338.56 

9.8053 

85000 

188.893 

4.4568E-001 

8 . 2195E-006 

275.52 

9.5496 

85050 

188.796 

4.4177E-001 

8 . 1516E-006 

275.45 

9 . 5494 

85100 

188.698 

4. 3790E-001 

8 . 0843E-006 

275.38 

9 . 5493 

85150 

188.601 

4. 3405E-001 

8.0174E-006 

275.31 

9 . 5491 

85200 

188.504 

4. 3024E-001 

7.9511E-006 

275.24 

9 . 5490 

85250 

188.406 

4. 2646E-001 

7 . 8853E-006 

275.16 

9.5488 

85300 

188.309 

4. 2271E-001 

7.8201E-006 

275.09 

9.5487 

85350 

188.212 

4. 1899E-001 

7 . 7553E-006 

275.02 

9.5485 

85400 

188.114 

4. 1531E-001 

7 . 6910E-006 

274.95 

9 . 5484 

85450 

188.017 

4. 1165E-001 

7 . 6273E-006 

274.88 

9.5482 

85500 

187.919 

4.0802E-001 

7 . 5640E-006 

274.81 

9.5481 
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